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j. ( Many messenger models with realistic gaugino masses are based on meta-stable vacua. In this work 

we study the thermal history of some of these models. Analyzing R-symmetric models, we point 
out that while some of the known messenger models clearly prefer the supersymmetric vacuum, 
Qh' there is a vast class of models where the answer depends on the initial conditions. Along with 

the vacuum at the origin, the high temperature thermal potential also possesses a local minimum 

(T) . far away from the origin. This vacuum has no analog at zero temperature. The first order phase 

>. 

transition from this vacuum into the supersymmetric vacuum is parametrically suppressed, and 

m ' 

the theory, starting from that vacuum, is likely to evolve to the desired gauge-mediation vacuum. 
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[ — ■ We also comment on the thermal evolution of models without R-symmetry. 
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I. INTRODUCTION 



It is well-known that low scale supersymmetry (SUSY) breaking mediated to the Standard 



Model through gauge interactions 



A Is* 14, [5] can successfully solve the supersymmetric 
flavor puzzle. Moreover, models of gauge mediation are often calculable. Thus gauge- 
mediation models should be considered as compelling candidates for the UV completion of 
the Minimal Supersymmetric Standard Model (MSSM). 

It has been suggested a long time ago (see e.g. f|) that the SUSY-breaking minimum can 
be metastable. It was further argued in [7(] that the gauge-mediated vacuum is likely to be 
metastable because of the messenger structure. For additional early examples of meta-stable 



11] that metastable 



supersymmetric models see e.g. 18, 9, 10] . Later it was shown in ISS 
SUSY-breaking minima are indeed generic and can be found, for example, in massive SQCD 
in the free magnetic phase. 

Some recent developments even more strongly advocate the idea of meta-stability of the 
gauge-mediated vacuum. It was proven in [12] that if gauginos acquire realistic soft masses 
within a calculable model, the gauge-mediation vacuum must be meta-stable within IR 
effective theory. This feature leads us to the question of which vacuum of the theory is 
cosmologically preferred. 

The issue of the thermal evolution of supersymmetric models with both supersymmetric 
and non-supersymmetric vacua was intensively studied in [3, 14, 111 16] in the context 



of ISS. All of these studies agreed that the ISS minimum is indeed cosmologically prefer- 
able. Moreover, as it was emphasized in 13) the theory is likely to end up in the non- 
supersymmetric vacuum even if it begins its thermal evolution from the supersymmetric 
vacuum. 

Nevertheless ISS-like meta-stability is not the type of meta-stability which is required in 
order to get a viable phenomenology. The ISS vacuum is stable within IR renormalizable ef- 
fective theory, and the supersymmetric vacuum becomes manifest only if non-renormalizable 
corrections are considered. Hence one cannot get leading order gaugino masses in ISS 
even if R-symmetry is maximally broken by some deformation. The absence of leading 
order gaugino masses was explicitly shown in the ISS-based models of direct gauge media- 
I20I ]. In order not to exacerbate the little hierarchy problem we will further 
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consider in our analysis only those theories which have leading order gaugino masses. 



2 



Thus we are forced to consider models of a rather different structure than ISS; dynam 



ica 
in 



examples of models of that type have been introduced e.g. in [21| and more recently 



221 ] . These examples, as well as other models one can potentially consider, are particular 



dynamical realizations of a much broader class of models of gauge mediation, so called extra- 



ordinary gauge mediation (EOGM) j23j. EOGM supplies the most generic ansatz for any 
calculable R-symmetric model with F-term SUSY breaking one can consider. 1 The thermal 
history of the EOGM models is the focus of our current work. 

R-symmetry is an important guideline in analyzing the models of EOGM. Spontaneously 
broken R-symmetry, even if it is an approximate accidental symmetry of the IR effective 



theory, is well motivated by the Nelson- Seiberg argument [25j| and because it ensures that 
leading order soft terms are CP-conserving. R-symmetry enables clear classification of the 
EOGM models, and endows these models with some useful generic properties, elaborated 



m 



231. 



In this work we systematically study the thermal history of EOGM models. We first 
concentrate on a special subset of EOGM, assuming the "minimal" UV completion. Namely 
we assume that the dynamics of the gauge mediation messengers is also responsible for SUSY 



breaking and R-symmetry breaking 26|]. The results which we obtain do not apply to the 
most general EOGM models. Nonetheless, as we show later, they turn out to be helpful in 
analyzing models, which are more complicated than just the "minimal completion" . 



Analyzing various EOGM models we stick to the definitions of 23(, distinguishing three 
different types of these models. Here we briefly describe these three types, while more precise 
definitions will be given in section HT1 

• Type I models. These are models with the full-rank messenger mass matrix. Those 
models are the most straightforward generalization of O'Raifeartaigh models. Type I 
models necessarily have vanishing leading order gaugino masses. 

• Type II models. These are models with full-rank coupling matrix between the messen- 
gers and the pseudomodulus. The simplest possible example of such models is minimal 
gauge mediation. 



1 As has been recently discussed in [241 ] . D-terms are only generated from F-tcrms in dynamical models 
of SUSY breaking. Thus, usually D-terms are parametrically suppressed and we will not consider this 
possibility in our current study. 
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• Type III models. Neither mass nor coupling matrices in these models are full-rank 
matrices. Those models share some similarities both with type I and type II models. 

We argue that type II models are cosmologically disfavored if the Universe is reheated 
above the SUSY-breaking scale. 2 Hence, in order to render these models viable one can 
obtain a new stringent bound on the reheat temperature. Interestingly, this bound goes pre- 
cisely in the same direction as bounds on the reheat temperature from gravitino cosmology, 



if the gravitino is more massive than a ke V 28|, [29 ] . 



Studying the type I models we point out that at high temperatures they might possess two 
different vacua. One of them is a global vacuum at the origin of the field space and the second 
vacuum is local and can potentially emerge very far away from the global vacuum. While 
the messengers in the local vacuum are at the origin of the field space, the pseudomodulus is 
stabilized far away from its origin. This second vacuum has no analog in the zero temperature 
field theory, and it disappears when the temperature drops below the SUSY-breaking scale. 
In spite of the fact that the free energy of this vacuum is larger than the free energy at 
the origin, the first order phase transition from this vacuum into the vacuum at the origin 
is strongly parametrically suppressed. 3 We also point out a necessary condition for this 
vacuum to exist. We show that this condition holds in type I models, that it does not hold 
in type II models and that it might hold in type III models. 

Turning to the type III models, we demonstrate that if the vacuum far away from the 
origin exists, the thermal history of the Universe strongly depends on the initial conditions. 
If it starts the evolution from the origin, it will inevitably slide into phenomenologically un- 
desirable minimum (or runaway). On the contrary, if the Universe is trapped after inflation 
in the vacuum far away from the origin, its transition rate into the vacuum at the origin is 
highly suppressed. When this minimum disappears the system slides to the gauge-mediation 
minimum. 

Our analysis assumes that the reheating temperature is higher than the messenger scale 
and that the Universe cools down adiabatically. We restrict ourselves in our analysis to 



2 Regarding some similar models, this observation has recently appeared in 2j|. Here we are studying these 
models in more generality. 

3 Some examples of the type I models, e.g. the minimal O'Raifeartaigh model, have been analyzed in [l4|. 
While we agree with their results regarding the minimum near the origin, we point out that one more 
thermal minimum might exist, significantly changing the analysis of the thermal history of such models. 
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calculable and renormalizable models. We assume that the Kahler potential is canonical 
at tree level. For the purposes of the current analysis we also assume that the Standard 
Model (SM) interactions are much weaker than the messenger model interactions; the SM 
interactions are completely neglected in our current study. It would be interesting to study 
the more general setup in the future. 

Our paper is organized as follows. In the second section we survey the thermal behavior 
of different types of "minimally completed" EOGM while also reviewing some facts about 
thermal field theory which are necessary for our analysis. In that section we show that a new 
vacuum far away from the origin may exist in type I models and explain intuitively why one 
should expect these minima to show up. In the third section we use these results in order to 
explain the behavior of some dynamical UV completions of EOGM found in the literature. 
We also comment on the behavior of some models which are not R-symmetric. Finally, in 
the fourth section we conclude. Some technical details are relegated to the appendix. 



II. SURVEYING REHEATED EOGM MODELS 



In this section we will survey all three types of R-symmetric EOGM models and analyze 
the thermal behavior of each. We will closely follow the notation of [23]. Consider the 
following superpotential 

W = (XijX + my) WW - fX = M lj ip i (p j - fX . (1) 

X may acquire both supersymmetric and SUSY-breaking VEVs X = (X) + 6 2 F. Note that 
the scalar component of X is a pseudomodulus, which is stabilized at loop level. 

Hereafter we will assume that the SUSY-breaking scale / is much lower than the messen- 
ger scale m. This assumption will further allow us to consider the couplings in the hidden 
sector to be of order 0(1) without sacrificing the zero-temperature meta-stability, and justify 
the formal assumption which we made A ^> gsM- 

At finite temperature T the thermal potential of the model is given by [30I ]: 

t 4 r°° ( I 2 1 "i 2 (0~ \ 

V(T) = V{T = 0) ± — / fa 2 ln l T e +... (2) 

= V(T = 0) + V^ loov {T = 0) + V th . 
Plus and minus signs stand for bosons and fermions respectively, and summation over all 
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degrees of freedom is implied. The ellipses stand for higher loop terms. Note that V(T = 0) 
denotes the tree-level zero-temperature potential. 

Since we are interested in understanding orders of magnitude of the phase transition 
temperatures, we will use the expansion of the expression ([2]) at temperatures much higher 
than the masses of the particles: 




The contribution of the particles much heaver than the temperature is Boltzmann suppressed 
and can be safely neglected. 



A. Type II models 

We first consider the models with det m = and det A ^ (the matrices m and A are 
defined in (JT])). These models are clearly phenomeno logically interesting since unlike type I 
models they provide viable gaugino masses. The reason for this is that leading-order gaugino 
masses are proportional to 

m 1 / 2 oc <9xhidet M . (4) 



On the other hand it was proven in [23| that if m is a full rank matrix, then R-symmetry 
necessarily renders the matrix Ai to be A-independent. Thus the models with detm = 
are of special phenomenological interest. 

Before discussing the thermal history of these models let us first briefly discuss their 
vacuum structure at zero temperature. Since gaugino masses are non-vanishing at leading 
order, the SUSY-breaking vacuum is necessary metastable. This feature is manifest: since 
det m = 0, we necessarily have messengers whose masses are given solely by AX. For 
sufficiently small values of X these messengers become tachyonic. Along the moduli space 
of X we have a stable messenger spectrum only for X > X m i n for some specific value of 
Xmin, which is of order of magnitude ^J-. Following the tachyonic direction the theory will 
either be stabilized at some supersymmetric minimum or will slide down to a runaway. 

Now let us turn to the thermal history of such models. To simplify the analysis we bring 
the matrix A to the diagonal form such that the messenger mass matrix is 

Mij = XiSijX + rriij . (5) 



6 



It is straightforward now to calculate the thermal potential of this theory at temperatures 
T ^> m, where by m we mean a rough messenger scale. The resulting expression is 

T 2 ( - \ T 2 

Vth = const + — f \M 2 \<P\ 2 + \M 2 W?\ + — Tr(A f A)|X| 2 (6) 

It is also clear in light of the results of appendix [A] that this thermal potential must be 
independent of rriij and /. The global minimum of the finite temperature potential is at 
the origin of the field space. One can also verify that no other minimum emerges far away 
from the origin, and at temperatures T ^> m the minimum at the origin of the field space 
is unique. 

Since the origin is not a stable point of the zero temperature effective potential, we 
should estimate the temperature at which a phase transition occurs. The messengers acquire 
tachyonic masses of order y/f, so one expects that the second order phase transition will 
occur at a comparable temperature. More precisely one can estimate for A < 1 using 
expression (jBJ) 




Now we analyze the behavior of the minimum at (fi = Cpi = with X getting a VEV at 
the messenger scale. Once the system is stabilized at this vacuum, SUSY is spontaneously 
broken and the MSSM fields get their soft masses. We will further call this minimum the 
"EOGM vacuum" . Let Xeogm denote the value of X at that minimum at zero temperature. 
When the temperature is well above the messenger scale, the effective potential for X scales 
as V ~ T 2 X 2 and the EOGM vacuum does not exist. Once the temperature drops below 
XX eogm all the messengers near the point X = X EOGM become heavy (since A is a full rank 
matrix) and the EOGM minimum emerges. Nonetheless the number of thermalized degrees 
of freedom in the EOGM minimum is smaller than at the origin. In general, one can notice 
that the free energy of the EOGM minimum is bigger than the free energy of the origin (or, 
alternatively, the supersymmetric vacuum after the second order phase transition occurs) at 
any temperature. 

This claim can be explained as follow. It is easy to see from ([3]) that a vacuum with 
a larger number of light degrees of freedom has lower free energy. When the EOGM vac- 
uum is formed, the only light supermultiplet in this vacuum is X. It becomes heavy at a 
temperature of order y^J- On the other hand, the vacuum at the origin has more light 
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degrees of freedom: both X and some messenger pairs are light (such messengers exist since 
rank A > rank m). Once we get to the temperature of order /, the free energy of the EOGM 
vacuum is already governed by the temperature independent term, f 2 , which arises since 
EOGM vacuum breaks SUSY. Hence the thermal potential at the EOGM vacuum below the 
temperature (j7j) scales as 

V EOG m(T 2 </) = / a -^T 4 + ...«/ a . (8) 

Clearly this value is bigger than the thermal potential along the path from the origin to 
the supersymmetric vacuum. The key point is that along that path the SUSY-breaking 
terms f 2 is absent and the value of the thermal potential is governed by the number of the 
massless degrees of freedom, which is necessarily smaller than (JHJ). Thus, type II minimally- 
completed EOGM models are always cosmologically disfavored unless the reheat temperature 
is significantly smaller than the messenger scale. Therefore in this context we find a new 
bound on reheat temperature. 



B. Type I models 

In this subsection we discuss the models with det m 7^ and det A = 0. Since these 
models cannot produce reasonable phenomenology, we will analyze this scenario as a toy 
model. The results will be useful in the analysis of type III models. 

Before analyzing genuine type I models, consider as an example the minimal 
O'Raifeartaigh model (let the fields ip, (p be gauge singlets in this example): 

W = XX(if 2 - f) + mtpip . (9) 

It is straightforward to calculate the thermal potential for small X in this model and verify 
that at sufficiently high temperatures all the fields can be stabilized at the origin. 
Let us now analyze this model at very large values of X, namely in the range 

X > T, X > m . (10) 

The key point for understanding this regime is noting that at zero temperature one finds 
very light particles in the O'Raifeartaigh models for X ^> m. Integrating out <p in for 
large X one finds the effective superpotential 
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Notice that the mass of (p is suppressed, rather than enhanced, by X. Since the zero- 
temperature potential for X vanishes at tree level, the non-thermal mass of the pseudomod- 
ulus is also strongly suppressed, such that both ip and X are thermalized in the regime (TTOT) . 
The mass of ip becomes exactly zero when X gets to infinity, so we should expect in this 
range thermal runaway for X. In order that this regime be reliable we need both inequali- 
ties in (11 01) to hold. If the temperature exceeds X, the field (p is also thermalized and the 
runaway behavior is lost. 

Now let us ask whether this runaway is stabilized. The thermal potential for X is 



m 4 T 2 



v « ~ m> ■ (12) 

But the zero temperature potential is a monotonically increasing function. The leading 
order term of the one-loop effective potential is 

Vi-ioo P In |X| 2 (13) 

which will inevitably stabilize the runaway. The mass matrix for the field p> is positive 
definite in the regime (fTUl) if m 2 > f. Thus this simple O'Raifeartaigh model, besides the 
vacuum at the origin, gets a local minimum of the thermal potential far away from the 
origin, X ^> T. More precisely the potential is balanced at 

x - ~ (s) TT ' (14) 

The free energy of this local vacuum is bigger than the free energy of the vacuum at the 
origin, so we are supposed to check if the thermal transition rate from this vacuum to the 
origin is suppressed. It turns out that the transition rate for the first order phase transition 
to the vacuum at the origin is parametrically suppressed since X* ^> T. We will justify this 
statement in the subsection III C\ after extending our analysis to broader class of models. 

Now we are ready to generalize this idea to the type I EOGM models. Consider for 
example a model with N messenger pairs and with the following superpotential: 

W = m(i)<pi<pi + \( Li )X<pi<pi +1 - fX (15) 

with i running over the values 1 . . . N. At large values of X all of the messengers excluding 
ipN and ipi are supposed to be integrated out. After integrating out the heavy messengers 
one gets the following superpotential: 

w = -fx+ R^> Tx»-^ ■ (16) 



Note that the power of X in this expression can be determined, without any calculation, 
from R-symmetry considerations. 

Now we check if the messenger mass-squared matrix is positive definite at zero tempera- 
ture. The mass-squared matrix for the scalars is given by two 2x2 diagonal blocks of the 
form 

(( det to \ 1 _|_ I AT i \ det to / \ 

In, a J ±t* i) n ^ xN (i7) 

_i_ / jy i \ det to / / det to \ 1 / 

Since we are discussing the case where X is large, this matrix is not necessarily positive 



definite. This is of course nothing but a manifestation of the claim of [23[ that the moduli 
space of the type I EOGM might become unstable when X exceeds some X max . A Nonetheless, 
for N = 2 (or, alternatively, multiple decoupled species of N = 2) the matrix ([17]) is always 
positive definite. 

Consider the thermal potential. Up to 0(1) numbers the leading terms of the thermal 
potential are 

/detm\ 2 T 2 /detm\ 2 T 2 ,, .„ . .„ s , , 

v,h ~ (tp:J + (tptJ Jxr [M + [M ) ■ (18) 

As expected, we get runaway behavior for X, and thermal masses for the light messengers 
which are highly suppressed compared to the naive expectation. The minimum always exists 
in models with two messenger pairs; in models with more messenger pairs it might exist for 
some definite range of temperatures. In these models one should check if the matrix (fTTj) is 
positive definite for a given choice of parameters. 

Now we are ready to estimate the position of the vacuum. We will continue doing analysis 
for generic N, but one should keep in mind that only for iV = 2 (or the set of decoupled 
models of this type) the existence of this vacuum is guaranteed. The one-loop effective 
potential far away from the origin is given by ffl3|) . 5 The minimum of the thermal potential 
scales as 

x ^((^y 1 ' 2 ^)^ . (19) 



Consider now that we start the thermal evolution in the vacuum far away from the origin 
and gradually lower the temperature. The vacuum will disappear either if X* ~ ^det m and 



4 For detailed discussion of this instability see 31 1 . 

5 Now this expression may be imprecise since we have different coefficients A^. Of course the dominant 



contribution here will come from the largest value of A. 
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integrating some of the messengers becomes unjustified, or if the messengers at that vacuum 
become heavier than the temperature, whatever happens earlier. The second condition is 
stronger; the vacuum disappears at temperature: 



namely right below SUSY-breaking scale. 

Note also that only for N = 2 this vacuum survives up arbitrarily high temperatures. If 
the number of messenger pairs is larger than 2, we get an upper bound on the temperature 
demanding that matrix (fTTj) is positive definite. This upper bound reads 



Evidently for sufficient separation between the messenger scale and the SUSY-breaking scale, 
the window defined by (120]) and (I2T]) can be very large. 

One can generalize these considerations to other messenger models, not necessarily of 
type I. Assume that we have several sets of messengers such that the matrix M. is block 
diagonal. We note that a necessary (but not sufficient!) condition in order to get this 
behavior is to have some subset of messengers which fulfill the condition 



If all of our messengers fulfill this condition, we get an EOGM model of type I, but in general 
type III models can also exhibit this behavior. 

C. Type III models 

Now we turn to the most phenomenologically interesting part of our analysis. These 
models have both det A = and det m = and share the features of both type I and 
type II models. In particular these models are well motivated since they lead to non- 
vanishing gaugino masses at leading order. The vacuum structure of these models at zero 
temperature is also non-trivial: they might develop instability both for very small values of 
the pseudomodulus (as in the case of type II models) and for very big values (as in the case 
of type I). 




(20) 




(21) 



rank m > rank A . 



(22) 
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The thermal histories of these models also share similar features with both type I and 
type II theories. It is clear that the thermal potential possesses a global minimum at the 
origin. If the Universe starts its evolution from this minimum, it will undergo a second order 
phase transition at the temperature (j7j), as was described in subsection III A[ and slide into 
the supersymmetric vacuum or runaway. 

Nonetheless, sliding into the supersymmetric direction is not the only possibility. Since 
these models might have some set of messengers fulfilling the condition (1221) . we can expect 
that certain models will also posses a local thermal vacuum far away from the origin. Here 
we show an explicit example of a model which indeed exhibits this behavior. Moreover, we 
show that starting the evolution from this minimum one ends up in the EOGM minimum, 
rather than in the supersymmetric one. 

Consider a model with following messenger mass matrix: 

/ X'X \ 
m XX 
M XX 
\ m J 

This model was analyzed in [23j as an example of the "minimal completion" of type III 
EOGM. It was shown that for sufficient separation of scales m and M the pseudomodulus 
can be stabilized between these two scales. 6 

Consider now this model thermalized when X is sufficiently far away from the origin. 
The first pair of messenger decouples, three other pairs maintain the condition (|22p . We 
can expect that the vacuum far away from the origin will show up. In order to know the 
exact scaling of this vacuum we just substitute iV = 3 in the formulas which we derived in 
the previous subsection. We find from ffl9|) 



M 



(23) 



X* 



/a A \" 1/2 m 2 MT 
V4W A 2 7 



1/2 



(24) 



As explained in the previous subsection, the vacuum far away from the origin for N > 2 
exists only at certain range of temperatures. In our case this range is given by 




1/2 



< T < 



4^ X 2 f 



(25) 



6 As explained in [22j |. the separation of masses is a necessary condition for the radiative stabilization of 
moduli in this kind of model. 
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Once the vacuum dissipates, the modulus X is at the messenger scale. Substituting the lower 
bound on the temperature from fl25|) into (1241) we explicitly verify that when the vacuum far 
away from the origin dissipates, the value of X* still respects an inequality X* > X EO gm- 

In order to understand what subsequently happens to the system which evolves through 
the vacuum far away from the origin, we analyze the thermal potential at the SUSY-breaking 
temperature. All the messengers are stable at the messenger scale, they have masses of order 
M and the potential rises in these directions, such that the pseudomodulus direction X is 
the only relevant direction for this discussion. Since at the temperature T ~ y/J none of the 
particles besides the multiplet X are thermalized in the vicinity of the EOGM vacuum, the 
shape of the thermal potential is the same as it is at zero temperature. Namely the potential 
monotonically rises at X > X EO gm approximately as V ~ logX + . . ., and decreases for 
smaller X towards the local SUSY-breaking EOGM minimum. Therefore we conclude that 
when the vacuum far away from the origin dissipates, the system undergoes second order 
phase transition to the EOGM minimum. On the other hand, the super symmetric minimum 
is located near the origin. Since the shape of the potential is similar to the shape of the zero 
temperature potential, the supersymmetric vacuum is separated from the EOGM minimum 
by a barrier with height of order y/f, and the distance between two these minima is of order 
M. Consequently, a second order phase transition from the vacuum far away from the origin 
into the supersymmetric vacuum is impossible. 

Now let us briefly discuss the hypothetical possibility of a first order phase transition 
from the vacuum far from the origin. The rate for this phase transition is given by 



In our case the potential rises from the origin as a quadratic function of X, up to the values 
X ~ T. At this value the potential for X starts falling as a negative power. The three 
dimensional bounce action can be approximated through the triangular approximation 32] : 



Since we obtained X* ^> T, this rate can be made arbitrarily small by appropriate choice of 
parameters for any temperature when the vacuum far away from the origin exists. Note that 
by the same reason the first order phase transition rate from the EOGM minimum at the 
temperature a/7 ( or lower) to the supersymmetric vacuum is also parametrically suppressed. 
In this last case one should of course substitute M instead of X* into (|2"7|) . 



(26) 




(27) 
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III. COMMENTS ON DYNAMICAL MODELS 



A. R-symmetric Dynamical and Retrofitted Models 

In the previous section we studied models of messengers with the "minimal" completion. 
It is plausible that at low energies this is a correct description, while all the small mass 
scales are just retrofitted as in [33(. In this case the description of the thermal behavior of 
these models is exhaustively explained in the previous section. Nonetheless, we will also be 
interested in applying our results to models with extra matter, which reproduce the EOGM 
ansatz dynamically. 

Consider first the type III models. Very few dynamical examples of this kind of model 
are known, one of them presented, e.g. in the model of "uplifted vacua in SQCD" j^. 7 
One can wonder if the mechanism of "rescuing" the meta-stable vacuum through evolution 
far away from the origin is still available in this case. The answer is negative. The crucial 
difference between the model presented in [22j and the "retrofitted" type III model is that in 
the retrofitted models the masses are "replaced" by the VEVs of some gaugino condensates. 
These condensates are usually formed at very high temperatures, and at low temperatures 
the masses behave as fundamental parameters. On the other hand, in dynamical models, 
these masses do not emerge before the thermal vacuum at the origin is destabilized and the 
relevant fields acquire VEVs. This usually happens at relatively low temperatures. Above 
that temperature the system evolves as if it did not have any masses at all, in other words 
similarly to the type II models. For example in the "uplifted vacua model" the masses 
do not emerge before the Universe is cooled down below the messenger scale, and hence 
the mechanism described in subsection III CI for rescuing the EOGM vacuum is inapplicable. 
Therefore, a stringent bound on the reheat temperature of these models applies. 

The next question we address is the behavior of (R-symmetric) models where messengers 
are added "by hand" and they are not part of the SUSY breaking sector. To illustrate this 
kind of model, consider an R-symmetric model of messengers and SM singlets which couple 
to the pseudomodulus. If the origin of some field is destabilized at sufficiently high temper- 
atures, they must have vanishing R-charge. 8 Hence, at sufficiently high temperatures X has 



7 An earlier example was given in 2l|. The thermal behavior of these two models does not differ. 



We prove this statement in the appendix. 
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a minimum at the origin independent of the details of the SM singlets to which it couples. 
If the theory starts from this vacuum, either by choice of initial conditions or because it is 
the only vacuum of the theory, it will inevitably cool down to the supersymmetric vacuum. 
The possibility of utilizing the rescuing mechanism through the vacuum far away from the 
origin strongly depends on the details of the theory. If the rank of the matrix m exceeds 
the rank of the matrix A for any subset of messengers or SM singlets, the vacuum far away 
from origin might emerge, as discussed in section [Til As previously noted, the masses should 
emerge at sufficiently high temperatures. 

34j , which is R-symmetric and includes 



Applying these arguments e.g. to the model 
both "hidden sector" fields and messenger fields (which behave as in type II), we conclude 
that the rescuing mechanism through the vacuum far away from the origin does not work. 
Even though the hidden sector fields perfectly fulfill the condition (1221) . the masses form 
only at the SUSY breaking scale, which is not sufficient to ensure the vacuum far away 
from the origin. We note, however, that the interactions between the messengers and the 
pseudomodulus in 34j can we extremely suppressed. While the model with the low scale 
M2 will favor the supersymmetric minimum, unless the reheat temperature is lower than 
the SUSY-breaking scale, the thermal history of the model with high M 2 will be largely 
governed by the SM interactions, which are not taken into account in our current study. 
If the SM D-terms indeed prevail, the final answer can change since the temperature for 
the second order phase transition into the EOGM vacuum can exceed the phase transition 
temperature in the supersymmetric direction. 9 



B. Comments on models without R-symmetry 



In this subsection we briefly comment on messenger models which lack R-symmetry. As 



36, 



37T| . 10 In these models the 



an example we will analyze here models presented in [35 . 
messenger sector was coupled to the ISS model, while the messenger masses were added "by 
hand". 11 These models exhibit a peculiar behavior: unlike R-symmetric models they might 



9 I am grateful to S. Abel, J. Jaeckel and V. Khoze f or p ointing mexmt the importance of the SM interactions. 

10 A similar model has been recently introduced in [38(. Since in [38( one should take A <C 1, the SM effects 
may become important. 

11 These masses were further "retrofitted" in [37f but it does not affect the thermal behavior of the model. 
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be stabilized at high temperatures at a point which is not a zero of the det(AX + m). 12 
Consider for simplicity the following superpotential: 

W = XX(fi<fi + X'Xcpjcpj - mtfiipi - fX . (28) 

The messengers tp, tp are charged under the SM model while 0, <p are uncharged. 13 It is easy 
to see that the thermal potential of the model (1281) up to an overall numerical factor is 

VthOcT^N^XXf + N^XX-ml 2 ) . (29) 

At sufficiently high temperature the pseudomodulus is stabilized at 

X = XN - m . (30) 
X'N^ + XN V 1 ; 

Both X = and X = m are zeros of det(AX — m), but the stabilization is at some 
intermediate value between them. 14 

Now we trace the evolution of this vacuum when the temperature drops down. Of course, 
we assume that the scale m is much larger than the SUSY breaking scale (otherwise even 
the zero-temperature meta-stability may be spoiled). At the scale m the vacuum at (130]) is 
destabilized and the second order phase transition occurs. The zero temperature masses of 
the fields tp, tp and <ft, (ft for X — X* are (XX* — m) and X'X* respectively. Those particles 
which are more massive will stop affecting the thermal potential at higher temperatures since 
their contributions will be Boltzmann suppressed, and the theory will slide into the direction 
where the lightest particles attract it. We can draw the conclusion that the models of ISS 
with messengers are not automatically safe; their thermal history strongly depends on the 
ratio jij^- If XN<p "wins" , the theory will be attracted to the direction of the supersymmetric 
minimum, but since these models have no restrictions on the hidden sector, this potential 
problem can circumvented by choosing a big enough magnetic group in the ISS sector. 

Now let us consider a more generic picture and try to draw some conclusions. A thermal 
potential of an arbitrary messenger model without R-symmetry is a sum of quadratic func- 
tions with the origins sitting at the zeros of det(AA + m). Thus the minimum of the potential 



12 Note, that in the ISS model (X) = 0. 

13 Those fields may be multiplets of some other symmetry group, but it is irrelevant for our discussion. 

14 As a consequence of appendix El an R-symmetric O'Raifeartaigh model cannot exhibit such behavior 
unless the R-charge of X vanishes. 
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is likely to be somewhere between these zeros. At zero temperature each of these zeros, if 
it is simple, is necessarily an attractive point of the potential. I.e. there is some vicinity 
of each of these points where the potential for X attracts it to the zero, independently of 
whether this point has tachyonic messengers 39]. 15 Namely each of these zeros can govern 
the thermal evolution, and an exact answer is highly model-dependent, as illustrated in the 
example of ( 1281) . 



IV. CONCLUSIONS AND OUTLOOK 



In this paper we studied the thermal behavior of various gauge mediated models, mostly 
with, but also without R-symmetries. We showed that at high temperatures all R-symmetric 
messenger models possess a minimum at the origin of field space. Starting from this min- 
imum the models generally evolve to the supersymmetric vacuum. In order to avoid this 
pitfall one should demand that the reheat temperature is lower than the messenger scale 
and that the Universe is trapped in the false vacuum after inflation. Demanding reheat 
temperature lower than the messenger scale can be quite a severe constraint on low-scale 
gauge-mediated models. 

We also showed that this stringent demand of the low reheat temperature can be circum- 
vented in certain models, where some particles, coupling to the pseudomodulus maintain 
the condition (1221) . In this case an additional minimum can emerge far away from the origin 
which can drive the thermal evolution of the Universe. In this case the Universe may end 
up in the SUSY-breaking minimum even for relatively high reheat temperature. Clearly the 
thermal evolution of models of this class strongly depend on the initial conditions after infla- 
tion. The subject of the initial conditions was not studied here and it would be interesting 
to study this subject in future works. 

We have presented convincing evidence that the thermal history is more special if the 
model is R-symmetric. It would be interesting to study further the connection between 
R-symmetry and the structure of the thermal potential. 

We also notice that we did not study the effects of the SM gauge interactions, concentrat- 
ing solely on the messenger sector analysis. While the messenger sector effects are extremely 

15 In a generic model without R-symmetry we expect all zeros to be simple. 
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important, and the thermal history of lots of models can be determined only by analyzing 
the messenger dynamics, it would be interesting to accomplish this analysis including the 
SM effects. Clearly the approximation where the SM effects are completely neglected cannot 
be valid in all models and in all parts of the parameter space, so it would be important to 
understand in the future, where precisely these effects are important and whether they can 
relax some of the constrains on the gauge-mediated models. 

Acknowledgements: I am indebted to Z. Chacko, Amit Giveon, Raman Sundrum, Brock 
Tweedie and especially Zohar Komargodski and David Shih for illuminating discussions. 
I am also grateful to Steven Abel, Joerg Jaeckel and Valya Khoze for the communication 
and discussion of the SM effects. The work was partially supported by NSF under grant 
PHY-0801323. 

APPENDIX A: PROPERTIES OF THE THERMAL POTENTIAL IN R- 
SYMMETRIC MODELS 

In this appendix we consider a thermalized R-symmetric O'Raifeartaigh model. We prove 
that at sufficiently high temperatures the origin of the field is destabilized only if its R-charge 
vanishes. We also show that the field can possess vanishing thermal mass at the origin only 
if there is more than one field in the theory with precisely the same R-charge. 

Consider the most generic renormalizable R-symmetric Wess-Zumino model: 

W = fa + ^-Mj + ^fo.op, . (Al) 

The couplings m^- and are symmetric in all indices. 

Now consider temperatures much higher than the masses in flAlj) and derive the lead- 
ing order field dependent thermal potential. Analyzing the approximation (J3]) we notice 
an important feature which significantly simplifies our analysis: the leading order thermal 
potential is merely sensitive to the trace of the mass squared matrix. Namely, we will not be 
be interested either in the SUSY-breaking terms in the zero temperature potential or in any 
other off-diagonal terms in the mass squared matrix. 16 In order to obtain a correct leading 

16 To use this property we need an additional assumption y/J -C m. Fortunately it holds in the vast majority 
of realistic gauge-mediated models. 
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term thermal potential it is sufficient to take a sum over the terms in the zero-temperature 
tree-level potential of the form 

V D ml^U* • (A2) 

The quadratic terms in the thermal potential are always positive since they are formed from 
the interaction in ( 1 All) . So the only terms which can cause destabilization in the thermal 
potential are either tadpoles or non-diagonal quadratic terms. 

In order to get a tadpole in the thermal potential we should have trilinear term in the 
zero temperature potential. The only way to get these interactions is 

V D m* kl 4>l Xijk4>i4>j + c.c. (A3) 

Since terms which may give contributions to the thermal potential should be of the 
form ( 1A2I) . either index i or index j should be equal to I. Namely the field <fij will have 
a tadpole in the thermal potential only if the terms 4>k<f>i and 4>k4>i4>j are both present in the 
superpotential. This can happen only if the R-charge of <pj is zero. 

Another source of thermal instability can potentially come from non-diagonal masses in 
the thermal potential. Such masses emerge from the quartic terms in the zero-temperature 
potential: 

VD A^A^I^IV^+c.c. (A4) 

But if k ^ I this necessarily means that the fields 4>k and <j>i possess the same R-charge. 
Working out the eigenvalues of the thermal mass matrix for <^ and (pi we conclude that one 
of the modes has positive thermal mass while the second is exactly massless. If this happens 
the fate of the massless mode depends on its zero-temperature mass. 
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